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Abstract 

Spin-polarized current effect is studied on the static and dynamic mag- 
netization of the antiferromagnet in a ferromagnet-antiferromagnet nano- 
junction. The macrospin approximation is generalized to antiferromag- 
nets. Canted antiferromagnetic configuration and resulting magnetic mo- 
ment are induced by an external magnetic field. The resonance frequency 
and damping are calculated, as well as the threshold current density cor- 
responding to instability appearance. A possibility is shown of generating 
low-damping magnetization oscillations in terahertz range. The fluctu- 
ation effect is discussed on the canted antiferromagnetic configuration. 
Numerical simulation is carried out of the magnetization dynamics of the 
antiferromagnetic layer in the nanojunction with spin-polarized current. 
Outside the instability range, the simulation results coincide completely 
with analytical calculations using linear approximation. In the instabil- 
ity range, undamped oscillations occur of the longitudinal and transverse 
magnetization components. 



1 Introduction 

The discovery of the spin transfer torque effect in ferromagnetic junctions 
under spin-polarized current [TJ [2] has stimulated a number of works in 
which such effects were observed as switching the junction magnetic con- 
figuration , spin wave generation '4 , current-driven motion of magnetic 
domain walls [5], modification of ferromagnetic resonance [6], etc. It is 
well known that the spin torque transfer from spin-polarized electrons 
to lattice leads to appearance of a negative damping. At some current 
density, this negative damping overcomes the positive (Gilbert) damping 
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with occurring instability of the original magnetic configuration. The cor- 
responding current density is high enough, of the order of 10^ K/cw?. 
This, naturally, stimulates attempts to lower this threshold. Various ways 
were proposed, such as using magnetic semiconductors [7|, in which the 
threshold current density can be lower down to 10^-10® A/cm^ because 
of their low saturation magnetization. However, using of such materials 
requires, as a rule, low temperatures because of low Curie temperature. 
Besides, the ferromagnetic resonance frequency is rather low in this case. 

In connection with these difficulties, the other approaches were pro- 
posed, based on high spin injection [8j or joint action of external magnetic 
field and spin-polarized current (5] [TD] . It seems promising, also, using 
magnetic junction of ferromagnet-antiferromagnet type, in which the fer- 
romagnet (FM) acts as an injector of spin-polarized electrons. The anti- 
ferromagnetic (AFM) layer, in which the magnetic sublattices are canted 
by external magnetic field, may have very low magnetization that pro- 
motes low threshold [11]. The AFM resonance frequency may be both 
low and high reaching 10^^ s~^, i.e. terahertz (THz) range. However, 
investigation and application of THz resonances is prevented because of 
their large damping. Such a damping in ferromagnetic junctions can be 
suppressed, as mentioned above, by means of spin-polarized current. The 
question arises about possibility of such a suppression in FM/AFM junc- 
tions. Note, that this problem has been paid attention of a number of 
authors [1^-^!]. 

Another interesting feature of the FM/AFM junctions with spin-polarized 
current is the possibility of canting the AFM structure by spin-polarized 
current without magnetic field. 



2 The equations of motion 

Let us consider a FM/AFM junction (Fig. [l| with current flowing perpen- 
dicular to layers, along x axis. An ultrathin spacer layer is placed between 
the FM and AFM layers to prevent direct exchange coupling between the 
layers. An external magnetic field is parallel to the FM magnetization 
and lies in the layer plane yz. The simplest AFM model is used with two 
equivalent sublattices. 

The AFM energy (per unit area), with uniform and nonuniform ex- 
change, anisotropy, external magnetic field, and the sd exchange inter- 
action of the conduction electrons with the magnetic lattice taking into 
account, takes the form [22] 

-/3'(Mi ■ n)(M2 • n) - ((Mi + M2) • H) 

-a»d((Mi+M2)-m) i, (1) 
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Figure 1: Scheme of the ferromagnet (FM)-antiferromagnet (AFM) junction; 
NM being a nonmagnetic layer. The main vector directions are shown. 

where Mi, M2 are the sublattice magnetization vectors, A is the uniform 
exchange constant, a, a' are the intra- and inter-sublattice nonuniform 
exchange constants, respectively, /3, /3' are the corresponding anisotropy 
constants, n is the unit vector along the anisotropy axis, H is the ex- 
ternal magnetic field, m is the conduction electron magnetization, Qsd is 
the dimensionless sd exchange interaction constant. We do not include 
demagnetization term because its contribution is small compared to the 
uniform exchange. The integral is taken over the AFM layer thickness 
Laf m ■ We are interested in the spin-polarized current effect on the AFM 
layer, so we consider a case of perfect FM injector with pinned lattice 
magnetization and without disturbance of the electron spin equilibrium, 
that allows to not include the FM layer energy in Eq. ((T|. 

Two mechanisms are known of the spin-polarized current effect on 
the magnetic lattice, namely, spin transfer torque (STT) [TJ [2] and an 
alternative mechanism |23l 124) due to the spin injection and appearance 
of nonequilibrium population of the spin subbands in the collector layer 
(this is the AFM layer, in our case). In the case of antiparallel relative 
orientation of the injector and collector magnetization vectors, such a state 
becomes energetically unfavorable at high enough current density, so that 
the antiparallel configuration switches to parallel one (such a process in 
FM junction is considered in detail in review [23). The latter mechanism 
is described with the sd exchange term in Eq. As to the former 

mechanism, it is of dissipative character (it leads to negative damping), 
so that it is taken into account by the boundary conditions (see below), 
not the Hamiltonian. 

The equations of the sublattice motion with damping taking into ac- 
count take the form 



dt 




dt 



+ 7 



[m. X H^^.] =0 (i = I, 2), (2) 
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where Mo is the sublattice magnetization, k is the damping constant, 



H 



a = 1, 2) 



(3) 



are the effective fields acting on the corresponding sublattices. 

From Eqs. ((T])-((3]) the equations are obtained for the total magnetiza- 
tion M = Ml + M2 and antiferromagnetism vector L = Mi — M2: 



dM 



1 K 

2 Mo 



M X 



dM 



+ 



T 



+7 [M X H] + 7 [M X H,d] 

+ ^7(/3 + I3'){M ■ n)[M x n] + ^-yiP - /3')(L • n)[L x n] 



+ -7(0 + a') 



M X 



d^M 



L X 



0, (4) 



dL 

dt 



1 K 

2 Mo 



L X 



dM 

dt 



+ 



+7 [L X H] + 7 [L X H,d] - 7A [L X M] 

+ ^7(/3 + /3')(M ■ n)[L X n] + ^^(fi - /3')(L ■ n)[M x n] 



where 



+ ^7(0 + a') 



H,d(x) = 



L X 



^^M 



dx^ 



M X 



dx^ 



(M(x') • m(x')) dx' 



0, (5) 



(6) 



SM{x) JO 

is the effective field due to sd exchange interaction. This field determines 
the spin injection contribution to the interaction of the conduction elec- 
trons with the antiferromagnet lattice. 

To find Tlsdix) field, the conduction electron magnetization m(a;) is 
to be calculated. The details of such calculations are presented in our 
preceding papers |26l [9]. Here we adduce the result for the case, where 
the antiferromagnet layer thickness Lafm is small compared to the spin 
diffusion length I with the current flow direction corresponding to the 
electron flux from FM to AFM: 

m= (m-l- Am)M, Am = -^^^^ f M(0) • Mf V (7) 

bLafm V / 

where m is the equilibrium (in absence of current) electron magnetization. 
Am is the nonequilibrium increment due to current, M = M/|M| is the 
unit vector along the AFM magnetization, Mf is the similar vector for 
FM, is the Bohr magneton, e is the electron charge, r is the electron 
spin relaxation time, j is the current density. 

It should have in mind in varying the integral ((6]), that the electron 
magnetization m depends on the vector M orientation relative to the FM 
magnetization vector Mf- From Eqs. ([6]) and (O we have 



Hsd = asdmM + a^d 

+asd 



eLAFM 



M(0) ■ Mf M 



Mf5{x - 0). 



(8) 
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By substitution ((8} into Q and ([5]), we obtain 



dM 



1 K 

2 Mo 



M X 

+7 [M x H] + jasd 



dM 

dt 



+ 



L X 



9L 

dt 



[m X Mf] 5{x - 0) 



+ + /3')(M ■ n)[M X n] + - /3')(L ■ n)[L x n] 



+ -7(q + a) 



M x 



9a; 2 



L X 



= 0, (9) 



9t 



1 K 

2 iWo 



L X 



dM 
'dT 



+ 



M X 



dL 

dt 



h7 [L X H] + -yasd^^^^ [l X Mf] <5(x - 0) 



-7 A 



asdin asdfJ-BTQj 



M 



eL/ 



[M 



(m(0) 



Mf I [L X M] 



+ ^7(/3 + /3')(M ■ n)[L x n] + ^jil3 - /3')(L ■ n)[M x n] 



+ ^7(q + a') 



L X 



d^M 



dx"^ 



+ 27(a-a') 



M X 



dx'^ 



0. (10) 



3 The boundary conditions 

The equations of motion ((9} and (|10p contain derivative over the space 
coordinate x. Therefore, boundary conditions at the AFM layer surfaces 
X = Q and x = Lafm are need to find solutions. The way of derivation 
was described in Ref. in detail. The conditions depend on the electron 
spin polarization and are determined by the continuity requirement of the 
spin currents at the interfaces. 

The terms with the space derivative in Eq. (|9} may be written in the 
form of a divergency: 



-7(0 + a) 



M X 



d^M 



1 , 
+ 27(" - 



L X 



d_ 

dx 

d^M 

dx 



-7(0? + q') 



M X 



dM 



+ -7(q-q') 



L X 



(11) 



The 3m vector is the lattice magnetization flux density. 

Let us integrate Eq. ([9| over x within narrow interval < x < e with 
subsequent passing to e — >■ +0 limit. Then only the mentioned terms with 
the space derivative and the singular term with delta function will con- 
tribute to the integral. As a result, we obtain an effective magnetization 
flux density with sd exchange contribution at the AFM boundary x — +0 
taking into account: 



J.//(+0) = Jm(+0) + 7a,,^^^^ [m(+0) X Mf] 



(12) 
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The magnetization flux density coming from the FM injector is 

J(-0) = ^jMf. (13) 



The component Jy = ( J(-0) • M(+0) ) M(+0) remains with the elec- 
trons, while the rest, 



Jx = J(-0) - J|| = [Mf - M(+0) (Mf ■ M(+0)) } 

[m(+0) X [m(+0) X Mf]] , (14) 



' eAP ■' 



is transferred to the AFM lattice owing to conservation of the magnetiza- 
tion fluxes [Dll]. 

By equating the magnetization fluxes (|12p and (I14p . we obtain 

^ _ MsQ • Tat Tat A'/rJl . MS^Q , 

) J 



eKP ■ 



[m X [m X Mf]] - Ta.d^^^j [m x Mf] , (15) 



all the M vectors being taken at 2: = 4-0. 

Since the AFM layer thickness is small compared to the spin diffusion 
length and the exchange length, we may use the macrospin approximation 
which was described in detail in Ref. [9]. In this approximation, the 
magnetization changes slowly within the layer thickness. This allows to 
write 

9Ja/ 3m{Lafm) -3m{+Q) ^ Jm(+0) ^^g^ 
dx Lafm Lafm 

because the magnetization flux is equal to zero at the interface between 
AFM and the nonmagnetic layer closing the electric circuit, J m (Lafm) = 
0. This allows to exclude the terms with space derivative from Eq. ((9]). 
In the rest terms, ]SA{x, t) and L(a;, t) quantities are replaced with their 
values at a; = 0. Then Eq. ((9| takes a more simple form: 



dM Ik 

IT ~ 2 lio 



M X 

dt 



+ 



^ dL 
L X — 

dt 



+ 7 [M X H] 



where 



+ + /3')(M ■ n)[M X n] + - /3')(L ■ n)[L x n] 

+K [m X [m X Mf]] + P [m x Mf] = 0, (17) 

ehAFMA'P ehAFM 

The term with delta function does not present here, since it is taken into 
account in the boundary conditions. 

Now we are to use again the macrospin approximation to exclude the 
space derivatives from Eq. (|10p . too. 

Owing to known relationships [22] between M and L vectors, namely, 
= 4Mo^ and (M ■ L) = 0, we have the following conditions: 

{^^ dM\ / aL\ „ / dM\ / ^ aL\ „ 
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By substituting Eqs. JTO} and JTT]) in JTg]) we find that conditions (fT9)) 
are fulfilled if the terms in pOI) 



-7(q + a') 



L X 



satisfy the following equations: 

(X ■ M) + K (^L ■ [m X [m X Mf 



M X 



PL 



X 



(20) 



(X.L) 



0. 



[m X Mf]) = 0, 
(21) 



Let us decompose the considered X vector on three mutually orthog- 
onal vectors: 

X = aM + 6L + C7 [L X M] . (22) 



The substitution ^ in ^ gives a = K (L ■ Mf j - P (JL x M] ■ MfJ , 
6 = 0. As to c coefficient, it is a current-induced correction to the co- 
efficient of 7[L X M] term in Eq. (|10[1 . i.e., a correction to the uniform 
exchange constant A. Let us estimate the correction. Multiplying (|22p 
scalarly by [L x M] with (|20p taking into account gives 



M2L2 



"L X Ml 



-7(a + Q) ILx +-7(a- 
_(a_Q)_ (L- — 



M X 



dx^ 



dx'^ 



It is seen that c ~ a/L\pj^.j, while A ^ ot/aF', where a is the lattice 
constant 22 . Since Lafm 3> a, the mentioned correction to A may be 
neglected. 

As a result, Eq. (|10|) takes the form 



dL 



1 K 

2 Mo 



L X 



dM 

dt 



+ 



M X 



dL 

dt 



-f 7 [L X H] - 7A [L X M] 

+ /^^(M ■ n)[L X n] + U{p - /3')(L ■ n)[M x n] 



[l X [p 



M X Mf +P L X Mf 



0. 



(24) 



Here, A constant contains also the equilibrium contribution of the con- 
duction electrons —ascirn/M . 

Equations ((TTJ and (1241) are the result of applying the macrospin con- 
cept to AFM. It is shown that such an approximation may be justified 
formally for AFM layer. Earlier, it was justified for EM layers [l] [2] and 
generalized [5] with spin injection taking into account. The macrospin 
approach corresponds well to experimental conditions and simplifies cal- 
culations substantially. The terms with K coefficient in Eqs. (|17|l . (|24|) 
describe effect of STT mechanism, while the terms with P coefficient take 
the spin injection effect into account. 
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4 The magnetization wave spectrum and 
damping 

We assume that the easy anisotropy axis hes in the plane of AFM layer 
and is directed along y axis, the FM magnetization vector is parallel to 
the positive direction of z axis, the external magnetic field is parallel to z 
axis too (see Fig. [T]) . 

We are interesting in behavior of small fluctuations around the steady 
state M = {0, 0, Mz}, L = {0, Ty, 0}, i. e. the small quantities M^, My, M, 

Mz — Mzi Lx, Ly = Liy i^y, 

Let us project Eqs. (|17ll . 1)24^ to the coordinate axes and take the 
terms up to the first order. The zero order terms are present only in the 
projection of Eq. (|24l) to x axis. They give 

P P 
Hz + - Hz + - 

Mz = 3_ ^ 7 



A+i(/3-/3') ^ 



Ly = ±\/4:M§-Ml^±2Mo. (25) 

Note that the spin-polarized current takes part in creating magnetic 
moment together with the external magnetic field due to the spin injection 
induced interaction of the electron spins with the lattice |23II24) . which P 
parameter in Eq. (|25p corresponds to. Such an interaction leads to appear- 
ance of an effective magnetic field parallel to the injector magnetization. 
As a result, a canted antiferromagnet configuration may be created with- 
out magnetic field. However, such a configuration corresponds to parallel 
orientation of FM and AFM layers, M||Mf. As is shown below, the insta- 
bility does not occur with this orientation, so that an external magnetic 
field is to be applied to reach instability. 

With Eq. (|25p taking into account, the equations for the first order 
quantities take the form 

dMx 1 «: r ^ dMy -p dLz 1 , „ , ^ 

-^■yiP + P'W^My - ^j{l3 - l3')LyLz + KMzMx = 0, (26) 
dMy Ik. — dM, 



Mz—^~{-lHz+P)Mx + KMzMy = 0, (27) 



dt 2 Mo dt 

dMz I K y dLx ,1 ,r, aisT r n ^oc^ 
— + 2 Mo ''^^ + 2^^^ " ^ ^^^^^ ^ °' ^''^ 

_ 1 « Uy^ ^M^<^k\^^HzI^Mz = 0, (29) 
dt 2Mo\ ^ dt dt j ' Mz 
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The set of equations (|26p - pip splits up to two mutually independent 
sets with respect to {Mx, My, L^) and {Lx, Ly, M^). They describe two 
independent spectral modes, one of them corresponds to precession of the 
AFM magnetization vector around the magnetic field, while another to pe- 
riodic changes of the vector length along the magnetic field. We begin with 
the spectrum and damping of the first mode. We consider monochromatic 
oscillation with oj angular frequency and put M^, My, Lz ~ exp(— itjf). 
Then we obtain from Eqs. (ge)), (g^), ([3l| 

{-zoj + KM,) Mx + l^-yH, + P- ^-yiP + /3')M. - 

i7(/^-/3')-i^}l.i^.=0, (32) 

{-zu; + KMz) My-UHz+P^]^ | = 0, (33) 



{-ioj + KMz) Lz + + -/?')- } LyMx = 0. (34) 

Note that aforementioned additivity (in the algebraic sense, the sign 
taking into account) of the external magnetic field and the injection-driven 
effective field takes place not only in the steady magnetization (|25p . but 
also in the oscillations of the magnetization and antiferromagnetism vec- 
tors, so that both fields appear in Eqs. (|32|l . (I33|) "on an equal footing". 

Usually, A 3> P, P' . With these inequalities and stationary solu- 
tion (|25p taking into account we find the dispersion relation for the mag- 
netization oscillation 

(1 -I- Hi^)uj'^ + 2iuu! — loq = 0, (35) 

where 

LOO = \l2-i^HAHE + {KM.Y + iiH, + (36) 



u = n'yHE + KM,, (37) 

He ~ AA/o is the exchange field. Ha — {13 — /3')Afo is the anisotropy 
field. Formulae (|36[) and (|37p (without current terms KM, and P) coin- 
cide with known ones [HIIT]. At He lO'^-lO'^ Oe, Ha ~ 10^ Oe we 
have oscillations in THz range, ujo ~ 10^^ s~^. In absence of current the 
damping is rather high: at k ~ 10~^ 
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Let us consider the contribution of spin-polarized current to the fre- 
quency and damping of AFM resonance. At first we consider STT mech- 
anism effect [Tl[^. According to JTSJ and ((23, 

KM. = -/^i. (39) 

At Hz < 0, that corresponds to direction of the magnetic field (and, 
therefore, the AFM magnetization) opposite to the FM magnetization, 
this quantity is negative. The total attenuation becomes negative also 
(an instability occurs), if 

.^ e.7Mo|g^WM ^^.^^ (40) 

At K ~ 10-^ -jMo ~ 10^° s-\ \Hz\ ~ 10^ Oe, Lafm ~ 10"'* cm, 
Q '-^ 1 we have jo ~ 10^ A/cm^. At j near to jo weakly damping THz 
oscillation can be obtained. At j > jo, instability occurs which may lead 
to either self-sustained oscillations, or a dynamic stationary state. The 
latter disappears with the current turning off. To answer the question 
about future of the instability it is necessary to go out the scope of the 
linear approximation. We have simulated numerically the behavior of 
the AFM magnetization behind the linear approximation (see section [8] 
below) . 

The spin-polarized current contributes also to the oscillation frequency. 
At the mentioned parameter values, we have \KMz\ ~ 10^^ s^^ that is 
comparable with the frequency in absence of the current. This allows 
tuning the frequency by the current or excite parametric resonance by 
means of the current modulation. 



5 Current-induced spin injection effect 

Now let us discuss the injection mechanism effect '23' '24'. As mentioned 
before, the role of the mechanism is reduced to addition of an effective 
field P/'y to the external magnetic field. At reasonable parameter values, 
that field is much less than the exchange field He, so that it does not 
influence directly the eigenfrequency (|36p . Nevertheless, that field can 
modify substantially the contribution of the STT mechanism, because 
Eq. (|39|) with 1)25^ taking into account now takes the form 

Such a modification leads to substantial consequences. At H. < 0, P < 
7|iJ2| the instability threshold H40|) is lowered, since \Hz\ — P/"f difference 
appears now instead of \Hz\. If, however, P > 7|-ffz| then the AFM 
magnetization steady state 

m ^ (42) 

becomes positive that corresponds to the parallel (stable) relative orienta- 
tion of the FM and AFM layers. In this case, the turning on current leads 
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to switching the antiparallel configuration (stated beforehand by means 
of an external magnetic field) to parallel one. With turning off current, 
the antiparallel configuration restores. 

Since the mentioned injection-driven field depends on the current (see (I18|l ). 
the instability condition (|4U|I is modified and takes the form 

<i<^, (43) 



l + rj rj 

where rj = a^d^^'lMoT, jo being defined with Eq. (|40[) . In absence of 
the injection mechanism, this condition reduces to H40[). Under rising 
role of this mechanism we have lowering the instability threshold, on the 
one hand, and the instability range narrowing, on the other hand. At 
j > jo/v the antiparallel configuration switches to parallel one. The 
relative contribution of the injection mechanism is determined with 77 
parameter. At typical values, agd ~ 10*, k ~ 10~^, 7M0 ~ 10^° s~^, 
T ~ 10~^^ s, this parameter is of the order of unity, so that the injection 
effect may lower noticeably the instability threshold. 

Now let us return to the set of equations (|26p - (|3ip and consider the 
second mode describing with Eqs. (|28[) - (|30p . The current influences this 
mode by changing steady magnetization Mz due to the injection effective 
field effect (see (|26p ). while the STT mechanism does not influence this 
mode. A calculation similar to previous one gives the former dispersion 
relation p5[) . but now 

u.-,^2j-H.H.-^, (44) 

.^.,H.:^. (45) 

At Hz < 0, P > \Hz I , that corresponds to current density j > jo /rj, the 
total attenuation becomes negative, while the frequency becomes imagi- 
nary, that means switching the antiparallel configuration to parallel one. 



6 Easy plane type antiferromagnet 

Let us consider briefiy the situation where AFM has easy-plane anisotropy. 
We take the AFM layer yz plane as the easy plane and x axis as the (hard) 
anisotropy axis. The magnetic field, as before, is directed along z axis. 

Without repeating calculations, similar to previous ones, we present 
the results. A formal difference appears only in Eq. (|36|) for the eigenfre- 
quency ujo of the first of the modes considered above. We have for that 
frequency 

ojo = \l{')Hz + PY + {KMzY. (46) 
The damping has the former form (|37p . so that the instability threshold 
is determined with former formula ((43]). 

In absence of the current [K = 0, P = 0) with not too small damping 
coefficient k, the frequency appears to be much less than damping, so 
that the corresponding oscillations are not observed. The current effect 
increases the frequency, on the one hand, and decreases the damping (at 
Hz < 0), on the other hand, that allows to observe oscillation regime. 
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7 Fluctuation effect 



It follows from Eq. (|43p that the threshold current density is proportional 
to the external magnetic field strength \Hz\ and decreases with the field 
lowering. A question arises about permissible lowest limit of the total field 
\Hz \ + -P/7- In accordance with Eq. (I25p . such a limit may be the field 
which create magnetization \Mz\ comparable with its equilibrium value 
due to thermal fluctuations. Let us estimate this magnetization and the 
corresponding field. 

The AFM energy change in V volume under canting the sublattice 
magnetization vectors with 6 < 180° angle between them is 

AE = AMo (1 - cos e)V = iAVMf , (47) 

the anisotropy energy being neglected compared to the exchange energy. 

The equilibrium value of the squared magnetization is calculated using 
the Gibbs distribution: 

(strictly speaking, the magnetization may be changed within (— 2Afo, 2Mo) 
interval, however, AVMq S> kT, so that the integration limits may be 
taken infinity). 

To observe the effects described above, the magnetization Mz which 
appears under joint action of the external field and the current (see (|25|) ) 
should exceed in magnitude the equilibrium magnetization (Mz)^^^. At 
the current density j = Jo/(l + ??) corresponding to the instability thresh- 
old, this condition is fulfilled at magnetic field 



\H.\> ^ — {l + v) = Hrmn. (49) 

At A ~ 10^, 77 ~ 1, Lafm ~ 10^® cm and lateral sizes of the switched 
element 10 x 10 /xm^ we have V ~ 10~^^ cm^ and Hmin ~ 30 Oe at room 
temperature. This limit can be decreased under larger element size. 

It should be mentioned also about other mechanisms of AFM canting. 
The most known and studied one is the relativistic Dzyaloshinskii-Moria 
effect (see, e.g. [22ii28j). Besides, possible mechanisms have been discussed 
due to competition between sd exchange and direct exchange interaction 
of the magnetic ions in the lattice [25]. At the same time, there are no 
indications, to our knowledge, about measurements of canting in conduc- 
tive AFM. So, present theory is related to conductive AFM, in which the 
lattice canting is determined with external magnetic field. 



8 Simulation 

With the purpose of simulating, it is convenient to modify slightly Eqs. p7|) 
and (|24p . namely: i) to use sublattice magnetizations Mi, M2 instead of 
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M, L, ii) to describe damping by the Landau-Lifsiiitz representation with 
double vector product, and iii) to introduce the following dimensionless 
variables: 

M,; = M,/Afo (j = 1, 2), h = H/Mo, r=-^^, 



1 + 

eLAFMfM^ 



With these variables, the set of equations (|17l) . (|24l) take the form 

- - [m, X h«^] ~ K [m, X [m, X h(;'J] , i = 1, 2, (50) 



dT 



h[\\. = h - AMa + /3 ( Ml ■ n ) n + /?' (Ma ■ n ) n + PqMp 



[fMi +M2) X Mf] 
-K"-^^^. —-2 ^=F(Mi, Ma), (51) 



Ml + M' 



(Ma, Ml) . (52) 



It is seen from Eqs. (|50|l - (I52|l that the motions of two sublattices are cou- 
pled each other. The sources of such a coupling are the uniform exchange, 
intersublattice anisotropy (coefficient /?'), as well the effective field due to 
the spin torque (the last term in Eq. ([ST])). 

We assume the following orders of magnitude of the used parame- 
ters: K ~ 10"^-10"^ A ~ 10^-10", /3, P' ~ 10"^ (/3 / 13'), Mo ~ 10^ 
G, 7M0 ~ lO^^s"^. Under such conditions, h — 1 corresponds to mag- 
netic field H ~ 10'^ Oe, and T — 1 corresponds to time t ~ 10"^" s. 
With asd ~ 10", T ~ 10"^^ s, Q ~ 1, Lafm ~ 10"'^ cm we have rjo ~ 10^ 
and Ko = 1 value corresponds to current density j ~ 10^ A/cm^. 

Equations (|50|) - (|52p written in coordinates represent six ordinary dif- 
ferential equations of the first order in the Cauchy form for Mix, Miy, 
Miz, M2x, M2y, M2z- The equations are not mutually independent be- 
cause of the normalization conditions 

|Mi| = |Ma|=l. (53) 

Nevertheless, all the six equations are used in our simulation, while the 
mentioned normalization conditions serve for checking correctness of the 
calculations. 

The simulation was carried out by means of Simulink program in 
MATLAB system with using Differential Equation Editor (DEE). Right- 
hand sides of the equations resolved with respect to derivatives were 
entered into the DEE block. The parameters A, /3, j3' , n, Kq, Pq were 
given as input signals, while three projections of the magnetization vec- 
tor Mx = Mix + M2x, My = Miy + May, M, = Mi, + M2, were 

output to oscilloscope blocks. Besides, jA/ij — (^Mix + Miy + Mi,^ 
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I " I 1-1 "2 " 2 \ ^/^ 

and M2 = I M^x + + M^z 1 values were output to digital dis- 
plays; these values must be equal to 1 (or, at least, be close to 1) under 
correct calculation. 

In the present work, we assume that the magnetic field and current are 
turned on at the initial time instant T = and hold constant. However, 
the procedure used allows to consider arbitrary time dependence of these 
quantities, specifically, to vary turning on and turning off instants. 

We began simulation with "verifying" results of the linear theory (ac- 
tually, this was a test for the model adequacy). As above, we assume 
that the AFM layer lies in yz plane, the current flows along x axis, the 
easy anisotropy axis coincides with y axis (n = {0, 1, 0}), the current is 
polarized along the positive direction of the z axis (Mf = {0, 0, 1}), the 
magnetic field is coUinear with z axis (h = {0, 0, ftz}). In such a configu- 
ration, the coUinear relative orientation of the FM and AFM layer magne- 
tizations is stationary (although, possibly, unstable), so that a small initial 
deviation from such orientation was given to imitate thermal fiuctuations. 
The initial value of the Mz component was chosen to be equal to the 
equilibrium value Mz = ftz/A in the given magnetic field without current. 
Thus, at A = 10^, /iz = -1 initial conditions Mi^ = M2z = -5 x 10"" 
(so that Mz = -1 X 10"^), Myx = = -5 x 10"*^, M^y = -M^y = 
(1 — M-xx — A^iz)^^^ are used taking the normalization conditions (|53p into 
account. 

With the dimensionless variables, the main results of the linear theory 
take the following form. 

The stationary magnetization in z direction is 

A> - M if fez + mKo _ fez + noKq 

Mz = Miz + A./2Z = ^^i(^,^,) - (54) 

Under deviation from the stationary value, two modes appear with di- 
mensionless (referred to 7M0) frequencies uii^i and damping i/i_2 defined 
with the following formulae: 



jI = 2(/3 - /3')A + (fez + ^oKof + ( , ) ' 



(55) 



.2-^ t;T^ (57) 

fez + ??oA^O 
ftz + 770^^0 

The instability of the antiparallel relative orientation of the FM and 
AFM layers at /iz < occurs in the current density value range defined 
with inequalities 

-^<i^o<^. (59) 
1 -I- KJ?o r\ii 

To start, the case of absence of the current [Kq = 0) has been con- 
sidered. Perfect agreement has been observed with Eqs. ((55)) - ((58)) . In 
particular, the magnetization oscillation frequency drops abruptly when 
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Figure 2: Oscillations of the longitudinal magnetization in the instability range 
at iCo = 7 X IQ-^. 

we put P — f}' , and the oscillations disappear completely, if we take, more- 
over, h = 0. At P ^ P' the observed frequency consists with Eq. (|55p . 
In absence of the magnetic field, the frequencies of the two modes coin- 
cide, so that the oscillation takes the form of a simple sinusoid. Under 
rising magnetic field, the frequencies become different, and beats appear 
because of interaction between the modes. The simulation results consist, 
also, with Eq. H54|) for the stationary magnetization in presence of the 
magnetic field. Turning on the magnetic field at T = instant leads to an 
aperiodic transient process which decays completely by T = 0.4, following 
which the magnetization remains constant value determined by Eq. (|54p . 

At A = 10^, K = 10"^, ryo = 10^, /iz = -1 the instability predicted by 
the linear theory must occur at Kq = 5 x 10~^ and disappear at Kq = 
1 X 10"^ In our numerical experiments, increasing Ko parameter from 
zero to the indicated threshold leads, in accordance with Eq. (|56p . to 
decrease of the damping of the magnetization vector precession about z 
axis because of the negative damping caused by the STT mechanism. 
Incidentally, the absolute value of the (negative) component decreases 
from 10~^ to 10~* because of infiuence of the injection mechanism, that 
consists, also, with Eq. (|54|) . Thus, the simulation results agree completely 
with the theory in the range below the instability threshold. 

Of course, the instability range (Kq > 5 x 10~^), which is not described 
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Figure 3: Beats of the transverse magnetization in the instability range at Kq = 
8 X 10-3. 

by the linear theory, is of more interest. At Kq = 5.1 x 10"'' undamped os- 
cillations are observed. At Kq — 7x 10~^ the precession oscillations of the 
transverse components Mx, My become almost sinusoidal with AT ~ 0.2 
period (this corresponds to angular frequency ~ 3 x 10^^ s~^ at the chosen 
parameter values) . The longitudinal component oscillates periodically 
with a negative stationary background, however, the oscillation form is far 
from sinusoidal one (Fig.[2|. At A'o = 8 x 10"^ the oscillations of M^, My 
components take on a form of beats (Fig. [3]), while they again become 
sinusoidal ai Ko = 9 x 10~^. At T^o = 1 x 10"^ (this is the right bound- 
ary of the instability range in the linear theory) all the three components 
oscillate around zero values. Further, at Kq ~ 1.1 x 10~^ the oscillations 
hold yet, but at Ko — 1.2 x 10"^ they disappear almost completely, and 
they are absent at Ko = 1.5 x 10~^. 

The subsequent increasing of the current leads only to rising magne- 
tization in the positive direction due to the injection mechanism. At Ko = 
2A/r;o, the longitudinal component Mz reaches the maximal possible value AIz = 
2 (the sublattices are flipped to a parallel position), then the (dimension- 
less) angular frequency is equal to 2A (this corresponds to ~ 10^^ s~^). 
Such a situation corresponds to rather high current density, it is estimated 
as 2 X 10* A/cm^ at the above-mentioned parameter values. 
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9 Conclusions 



The obtained results show a prineipal possibility of eontroUing frequeney 
and damping of AFM resonance in FM/AFM junctions by means of spin- 
polarized current. Under low AFM magnetization induced by an external 
magnetic field perpendicular to the antiferromagnetism vector, the thresh- 
old current density corresponding to occurring instability is less substan- 
tially than in the FM-FM case. Near the threshold, the AFM resonance 
frequency increases, while damping decreases, that opens a possibility of 
generating oscillations in THz range. 

Numerical simulation allows to trace behavior of the FM/AFM junc- 
tion in the whole current density range. The instability range predicted 
by the linearized theory is broadened only slightly because of nonlinear 
effects. In the instability range undamped oscillation of sinusoidal or more 
complicated form including beats. 

Under magnetic fields low compared to the exchange field, the induced 
magnetization is small in comparison with the sublattice magnetization, so 
that the stationary oscillation amplitude beyond the instability threshold 
is low, too. 

Thus, the following features may be expected in comparison with the 
similar effects in FM/FM junctions. First, the instability threshold is to 
be lower because of the lower magnetization. This is a favorable fact facili- 
tating observations. Second, the oscillation intensity beyond the threshold 
also lowers as a square of the magnetization. This may make the effect 
difficult to observe. Nevertheless, studying the current-driven nonlinear 
oscillations in FM/AFM structure is of principal interest, because the cur- 
rent induced instability can occur at relatively low current density, ^ 10^ 
A/cm^ 

The simulation results reveal an interesting possibility of a spin-flip 
transition without magnetic field under the action of a high-density spin- 
polarized current only. Such a current overcomes the exchange forces and 
aligns the sublattice moments in parallel. Under such conditions, applying 
a low alternating magnetic field can excite precession of the magnetization 
vector at the AFM resonance frequencies, which may be as high as 3 x 10^^ 
s~^ or more. Such a TKIz resonator might be useful to detect and measure 
signals in THz range. 
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